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Abstract 

We report a search for the h c meson via the decay chain — > h c K^, h c — > r/ c 7 with rj c — ► 
KgK ± ir T and pp. No significant signals are observed. We obtain upper limits on the branching 
fractions for S ± — ► r] c ^K ± in bins of the rj c j invariant mass. The results are based on an analysis 
of 253 fb _1 of data collected by the Belle detector at the KEKB e + e~ collider. 

PACS numbers: 13.25.Hw, 13.20.-v 
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The h c meson is the l 1 Pi spin singlet state of cc, which is one of the missing states in the 
charmonium spectrum below DD threshold. It is expected to be a narrow resonance (I\ c < 
1.1 MeV/c 2 ) that lies between J/ip(lS) and ip(2S). The predicted masses of h c vary and are 
summarized in Ref . Q • A typical value is less than 10 MeV from the center of gravity of the 
l 3 Pj states (xo, X d and Xc2 ), which is M c . . g . = (M Xc0 + 3M Xcl + 5M Xc2 )/9 = 3525.4 ± 0.1 
MeV/c 2 . The h c meson should decay dominantly to r]^ with a branching fraction of about 
50% Qii- 

The E760 collaboration has reported an enhancement in the pp —>■ h c — > J/ipir cross 
section and identified it as the l x Pi state with a mass of 3526.2±0.25 MeV/c 2 @. This result 
was not confirmed by the subsequent experiment E835 with significantly higher statistics. 
However, E8350 reported promising evidence for the h c in h c — > rj c ^. Recently, CLEO0 
has reported the observation of h c — ► rjc'j at a mass of M = 3524.4 ± 0.6 ± 0.4 MeV/c 2 . The 
masses obtained by CLEO and E835 are within 1 MeV of M c o g . 

M. Suzuki and others have proposed using the decay chain B — > h r K, h c —>■ tvy 
to look for the h c meson. Other charmonium candidates including the ?7c(25')j8j, X(3872)[9j 
and y (3940) [13 were first observed in two-body B decays, where the kinematic constraints 
from the exclusive B decay and production at threshold provide substantial background 
reduction. The decay amplitudes for B — »• h c K and B — >• Xco,2K vanish in the factorization 
limit. The branching fractions for B + — > x c o,2-^ + have been measured and the results 
are given below [l3|, 



B(B + -> XcoK + ) = (1.34 ± 0.45 ± 0.15 ± 0.04) x 10" 4 (BaBar), 

B(B + ^XcoK + ) = (1.12 ± 0.12 ± 0.18 ± 0.08) x 10- 4 (Belle), 

B(B + -> Xc2K + ) < 0.3 x 10" 4 . (1) 

The fairly large branching fraction for B + — > Xd)K + suggests that non-factorizable contri- 
butions in B decays to charmonium can be sizable. The decay B — > h c K may occur via 
the color octet mechanism or re-scattering processes [16j at a rate comparable to that 
of the factorization allowed decay mode B — > XaK. Thus, measurement of the branching 
fraction for B — ► h c K will provide useful information on non-factorizable contributions in 
B to charmonium decays. 

Here we present the results of a search for B + — > h c K + , h c —>■ rjcj with % — > KgK ± 7T T 
and using a 253 fb _1 data sample, which contains 275 x 10 6 produced BB pairs. The 
data were collected at the T(4S) resonance with the Belle detector at the KEKB e + e~ 



collider [17J. In addition, we use a 28 fbr 1 data sample collected at an energy 60 MeV below 
resonance to measure the continuum background. 

The Belle detector is a large-solid-angle magnetic spectrometer that consists of a silicon 
vertex detector (SVD), a 50-layer central drift chamber (CDC), an array of aerogel threshold 
Cerenkov counters (ACC), a barrel-like arrangement of time-of-flight scintillation counters 
(TOF), and an electromagnetic calorimeter (ECL) comprised of CsI(Tl) crystals located 
inside a super-conducting solenoid coil that provides a 1.5 T magnetic field. An iron flux- 
return located outside of the coil is instrumented to detect K® mesons and to identify muons 
(KLM). The detector is described in detail elsewhere jT^. A sample of 152 million BB pairs 
was taken with a 2.0 cm radius beampipe and a 3-layer silicon vertex detector; another 
sample of 123 million BB pairs was taken with a 1.5 cm radius beampipe, a 4-layer silicon 
detector and a small-cell inner drift chamber 

Event selection criteria were determined using the figure of merit, which is defined as 
S/ \/~S + B, where B is the number of background events and S is the number of signal events 
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in a GEANT-based Monte Carlo simulation. We assume B(B + — > h c K + ) = 3 x 10 -4 and 
B[h c — > 7^7) = 0.5 for the signal and determine the background S from the sideband data. 
Signal events are simulated for five different values of the h c mass, which are M^ c = 3.23, 
3.33, 3.43, 3.527 and 3.63 GeV/c 2 , and assuming the intrinsic width Th c = 1 MeV/c 2 . We 
determine the final optimization of selection requirements with the M/ lc = 3.527 GeV/c 2 
MC sample. 

We select well measured charged tracks with impact parameters with respect to the 
interaction point (IP) of less than 0.3 cm in the radial direction and less than 5 cm in the z 
direction, which is opposite to the positron beam direction. The tracks are required to have 
the transverse momenta greater than 50 MeV/c and have more than 6 axial and 2 stereo 
CDC hits. 

Particle identification likelihoods for the pion and kaon particle hypotheses are calculated 
by combining information from the TOF and ACC systems with dE/dx measurements in 
the CDC. To identify kaons, we require the kaon likelihood ratio, Lk/{Lk + L n ), to be 
greater than 0.6, which is 89% efficient for kaons with a 8% misidentification rate for pions. 
For the charged kaons that come directly from the B meson rather than from the subsequent 
decay of the rj c , the kaon likelihood ratio is required to be greater than 0.5. In addition, we 
remove all kaon candidates that are consistent with being either protons or electrons. To 
identify pions, we require Lk/{Lk + L w ) to be smaller than 0.7, which is 94% efficient for 
pions with a 12% misidentification rate for kaons. 

Protons and antiprotons are identified using all particle identification systems and are 
required to have proton likelihood ratios [L P /(L P + L K ) and L p /(L p + L w )] greater than 
0.5. Proton candidates that are electron-like according to the information from the ECL are 
vetoed. This selection is 90% efficient for protons with a 6% misidentification rate for kaons 
and a 3% misidentification rate for pions. 

We select K s — ► 7r + 7r~ candidates from pairs of oppositely charged tracks that are consis- 
tent with the pion hypothesis to form common vertices and lie within the mass window 0.482 
GeV/c 2 < M(7r + 7r~) < 0.514 GeV/c 2 , which corresponds to ±4cr. The K s vertex is required 
to be displaced from the IP; the vertex direction from the IP is required to be consistent 
with the K s flight direction. The K s requirements are described in detail elsewhere jiol ]. 

We reconstruct r\ c candidates in the K s K ± n T and pp decay modes. The r\ c candidate is 
required to have an invariant mass in the range between 2.935 and 3.035 GeV/c 2 . In order 
to reduce the combinatorial background, the charged daughters of the r] c are required to 
come from a common vertex that is consistent with the interaction point profile. 

Photon candidates for the decay h c — > rjcj are selected from ECL clusters that are not 
associated with charged tracks extrapolated from the CDC. We require the photons have 
energies above 60 MeV, and at least five crystal hits. 

To isolate the signal, we form the beam constrained mass M bc = y -Ebeam — -Brecon anc ^ 
energy difference AE = E rccon - E heam , where -E bea m, -Brecon and P reC on are the beam energy, 
the reconstructed energy and the reconstructed momentum of a B candidate in the T(4S) 
center of mass frame. The signal region for M bc is 5.270 < M bc < 5.290 GeV/c 2 . The 
signal region for AE is —50 < AE < 35 MeV, which corresponds to ±2.5cx where a is 
the resolution determined from a Gaussian fit to the Monte Carlo simulation. If more than 
one signal candidate is found in an event, we select the one with the largest invariant mass 
M(K +r y) 1 where the kaon comes from the B, the best \ 2 °f the r] c vertex and the invariant 
mass M(i^°fT ± 7r T ) or M(pp) closest to the nominal mass of the rj c . These requirements are 
imposed in the order listed till only one candidate is selected. 
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To suppress the large background from continuum e + e" — > qq where q = u,d, s, c, we 
first remove events with the normalized second Fox- Wolfram moment R2 > 0.5. We then 
form a likelihood ratio using two variables. Six modified Fox- Wolfram moments jijj and the 
cosine of the thrust angle are combined into a Fisher discriminant T . For signal Monte Carlo 
and continuum data, we form probability density functions for this Fisher discriminant, and 
the cosine of the B decay angle with respect to the z axis {cos6 B ). We then calculate the 
likelihood ratio 1Z = Cs/{Cs + £bg) f° r the B + — > h c K + signal Monte Carlo and continuum 
data. The likelihood ratio TZ for the r] c —>■ KgK~n + mode is required to be greater than 
0.7. This requirement retains 70% of the signal while removing 92% of the continuum 
background. For the pp mode, which has less continuum background, we require 1Z to be 
greater than 0.6. 

In addition to backgrounds from continuum there are also backgrounds from other B 
decays. To investigate these backgrounds, we use a sample of 379 x 10 6 BB Monte Carlo 
events. We find that the dominant backgrounds come from B + — > t] c K* + , K* + — > K + tt° 
and B° — > rj c K*°, K*° —> K + tc~. These backgrounds peak in the Mb c distributions. We 
reject the events if the photon combined with any other photon makes a tt° candidate with 
0.114 < M(77) < 0.151 MeV/c 2 . This tt° veto requirement is 83% efficient for signal and 
removes 51% of background 7r°'s. We also require the cosine of the angle in the ^7 rest 
frame between the 7 and the kaon coming from the B candidate be smaller than 0.6 (0.9) 
if the invariant mass M{ri c r )) is smaller (greater) than 3.5 GeV. These requirements retain 
70% (85%) of the signal while removing 67% (72%) of the B — > r] c K* backgrounds. 

We also find backgrounds from B + — > i] c K + and B — > _£)(*)+£)(*) ; which peak in the 
Mb c distributions. We remove the B + — > r] c K + background if the AE for this decay mode 
is between —60 and +60 MeV (±6cr). This requirement is about 100% efficient for the 
signal. We also apply a Df veto if the invariant mass M(K + Ks) is in the range 1.938 < 
M(K + K S ) < 1.998 GeV/c 2 (±3<r). This requirement is 94% efficient for the signal. _ 

After all the event selection requirements are applied, no significant peaking BB back- 
grounds are observed. A component for the BB background is included in the fits to data. 
In addition, uncertainties in the BB background composition are included in the systematic 
error. 

In Fig. H we show the M^ c and AE distributions in five 100 MeV/c 2 bins of M(t7 c 7), 
which correspond to a ±3a range around the central values used in the signal MC. We 
determine signal yields from unbinned two-dimensional maximum-likelihood fits to the M^ c - 
AE distributions in the in region 5.2 < M bc < 5.3 GeV/c 2 and -0.2 < AE < 0.2 GeV. 
We use a Gaussian function for M^ c and a double Gaussian for AE to model the signal. 
The mean of the Gaussian for Mb c is fixed to 5.279 GeV/c 2 , which is determined from a 
B + — > D*°7T + data sample, and other parameters are fixed to the values from MC simulation. 
The functions used to model the backgrounds for other B decays are determined from MC. 
The continuum background is modeled with an ARGUS background function that behaves 
like phase space near the kinematic boundary [22] for M^ c and a linear function for AE. We 
determine the shape parameters for the background functions from a fit to the off-resonance 
data sample. We find that the background shapes do not depend on M{rj c ^). Therefore, we 
combine the data in the range 3.17 < M{rj c ^) < 3.72 GeV/c 2 to increase statistics. 

Because the mass of the h c is not well established, we fit the M^ C -AE distributions in the 
five 100 MeV/c 2 bins of M(t] c j). The results of the fits are shown in Fig. ^ The signal yields 
and the detection efficiencies, which are determined from the signal MC samples described 
above, are given in Table HI No significant signals are observed for 3.17 < M(j] c j) < 3.67 
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FIG. 1: The Mt, c distributions in the AE signal region and the AE distributions in the M^ c signal 
region for the rj c — > K^K'ir^ (left) and r\ c — > pp modes (right) in 100 MeV bins of M(rj c j) for 
3.17 < M(r/ C 7) < 3.67 GeV/c 2 . The distributions are shown in the increasing order of M(r/ C j) from 
the top to the bottom. The solid curves are the results of the fits. The dashed curves represent 
background components from B decays. 



To check for possible binning effects, we determine the branching fractions for the 7/ c 7 
invariant mass ranges that are shifted by 50 MeV/c 2 with respect to the nominal range. The 
results of these fits are shown in Fig. |21 The signal yields and the detection efficiencies are 
given in Table |l] No significant signals are observed for the range 3.22 < M(^ c 7) < 3.72 
GeV/c 2 . 

Fig. El shows the M{j] c ^) distribution for data in the Mb c and AE signal region (points 
with error bars). The distribution is consistent with the background determined from the 
Mb c sideband data in the region 5.20 < Mb c < 5.26 GeV/c 2 . The expected contribution for 
a resonance of a mass of 3.527 GeV/c 2 with a branching fraction at the observed upper limit 
is also shown. 

To demonstrate the effectiveness of the analysis procedure and the method for branching 
fraction determination, we examine the decay chain B + — > xaK + , Xci ~^ J/'i'li J/^ ~ > PP- 
A clear signal of 14. 8^3^ events is observed in the M^-AE distribution. We use the yield and 
the MC detection efficiency of 0.131 to determine the branching fraction B{B + — > xaK + ) = 
(6.1l]|) x 10" 4 , where the error is statistical only. This is in very good agreement with the 
result calculated from the world averages |l2j, B(B + — > Xc\K + ) = (6.8 ± 1.2) x 10~ 4 . 

The largest contributions to the systematic error in the detection efficiency are the un- 
certainties in the efficiencies for tracking, particle identification and photon detection. The 
errors in kaon, pion and proton identification efficiencies are obtained from kinematically 
selected D* + — > D°7T + , D° — > K~n + and A — > pn~ decays in the data. We apply correction 
factors of 0.975 and 0.941 for the pion and proton detection efficiencies, respectively. The 
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FIG. 2: The Mb c distributions in the AE 1 signal region and the Ai? distributions in the M\, c signal 
region for the rj c — > KgK~TT + (left) and ?] c — > pp modes (right) in 100 MeV bins of M(j] c ^) for 
3.22 < M(r/ C 7) < 3.72 GeV/c 2 . The distributions are shown in the increasing order of M(r/ C j) from 
the top to the bottom. The solid curves are the results of the fits. The dashed curves represent 
background components from B decays. 



TABLE I: Detection efficiencies and signal yields in 100 MeV/c 2 bins of M(r/ C 7). 

rj c -> K^K~ir + r] c -> pp 

Mjrjc-y) (GeV/c 2 ) e (%) Yield e (%) Yield 

3.17- 3.27 4.39 ±0.10 — 3.2±|;| 11.7 ±0.2 = 4Ji|f 

3.27- 3.37 4.62 ±0.10 7.0±f;§ 12.6 ± 0.2 — 1.2±^| 

3.37 - 3.47 5.31 ± 0.10 — 3.6±|;| 15.1 ± 0.2 -0.30t 2 ;? 

3.47 - 3.57 5.47 ± 0.10 — 3.9±|;| 15.1 ± 0.2 -0.761};? 

3.57- 3.67 5.95 ±0.11 3.2±|;g 16.6 ± 0.2 1.4l 2 ;° 

3.22 - 3.32 4.66 ±0.10 -2.2l"£o 12.2 ±0.2 -3.8t 2 ;g 

3.32 - 3.42 4.83 ±0.10 0.70±|;5 13.8 ± 0.2 1.4±f;| 

3.42 - 3.52 5.41 ±0.10 2.1±|;| 14.9 ± 0.2 1.30±f*| 

3.52 - 3.62 5.48 ±0.11 — 1.0±|;| 16.2 ±0.2 OT^ 5 

3.62 - 3.72 6.02 ±0.11 0.09±|;g 16.9 ± 0.2 l.Ol 2 ;^ 



particle identification systematic error is 5.1% for the r] c — > KgK~7r + mode and 6.6% for the 
Vc PP mode. The uncertainty in photon detection efficiency is determined to be 5% using 
a V ~~ * 77 data sample. The uncertainty in the r\ c vertex reconstruction is estimated to be 
2% using a <p ~~ > K + K~ sample. The systematic error due to the modeling of the likelihood 
ratio cut is determined to be 4% using B + — > D n + events reconstructed in data. We also 
include the MC statistical uncertainty and the uncertainty in the number of BB pairs in the 
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FIG. 3: The M(K^K ± 7t^'y) distributions. The points with error bars are the data in the My, c and 
AE signal region; the solid histogram is the M\, c sideband data in the region 5.20 < Mb c < 5.26 
GeV/c 2 ; the hatched histogram shows the expected contribution for a signal with a branching 
fraction at the observed upper limit for a resonance with M(rj c j) = 3.527 GeV 

data sample. The sources of systematic error are combined in quadrature to obtain the final 
systematic error in the detection efficiency, which is 10.3% for the r] c — > KgK~ir + mode and 
10.1% for the r\ c — > pp mode. 

The uncertainty in the signal yield from the fit is determined by varying the mean of 
the signal for M^ c by 0.5 MeV/c 2 , and all other shape parameters of the signal and the 
background by la of the measured errors. The results are combined in quadrature to obtain 
the total uncertainty, which depends on M(^ c 7) bin and ranges from ±0.5 to ±2.1 for the 
1] C — > KgK~ir + mode and from ±0.0 to ±1.3 for the rj c —>■ pp mode. 

We combine the likelihoods for the r] c — > KgK~ir + and rj c —> pp modes, taking into 
account the respective systematic errors in the detection efficiencies and uncertainties in the 
signal yields. We determine upper limits at 90% confidence level (C.L.) on the branching 
fractions for B + — ► r) c r yK + . The results are given in Table [H] in bins of M{j] c ^). 

In summary, we have searched for the h c meson in the decay chain B + — > h c K + , h c — > rjcj, 
where the r\ c is reconstructed in the KgK ± ir T and pp modes. No significant signals are seen 
for 3.17 < M( / ~ff] c ) < 3.67 GeV/c 2 . We obtain upper limits on the branching fractions for 
B + — > ^yi] c K + for different ^7 invariant mass ranges. Assuming B(h c —>■ / ~fr] c ) = 0.5, these 
results give 90% C.L. upper limits on branching fractions for B + — > h c K + as a function of 
the h c mass. The results are shown in Fig. For Mh c = 3.527 GeV/c 2 , we find B{B + — > 
h c K + ) < 3.8 x 10~ 5 . This is below the lower bound on the B — > h c K branching fraction 
obtained by Colangelo, Fazio and Pham [la ], which is B(B — > h c K) = (2 — 12) x 10 -4 . These 
results are comparable to the upper limit for B + -> x c2 if + [Ll 

but below the measured rate 

for B + — > XcqK + , two other non-factorizable decays. The upper limits obtained in this 



TABLE II: Upper limits at 90% C.L. on branching fractions for B + — > ^yrj c K + in bins of M(r] c ^f). 

M{rjc~)) (GeV/c 2 ) Branching Fraction 



3.17- 


3.27 


< 5.9 x 10" 


-5 


3.27- 


3.37 


< 8.6 x 10" 


-5 


3.37- 


3.47 


< 3.2 x 10" 


-5 


3.47- 


3.57 


< 3.8 x 10- 


-5 


3.57- 


3.67 


< 5.8 x 10" 


-5 


3.22 - 


3.32 


< 4.7 x 10" 


-5 


3.32 - 


3.42 


< 6.7 x 10" 


-5 


3.42 - 


3.52 


< 6.2 x 10" 


-5 


3.52 - 


3.62 


< 2.8 x 10" 


-5 


3.62 - 


3.72 


< 3.9 x 10" 


-5 



paper assume B(h c — > 7%) = 0.5 and therefore must be renormalized when this h c absolute 
branching fraction is measured. These results may also be used to constrain branching 
fractions of other charmonium or charmonium-like states that decay to r\(Pj. 
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